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ABSTRACT 


The response of a closed cylindrical shell is determined for acoustic 
excitation which is random in time but deterministic in space. Two slightly 
different formulations of the statistical energy method are utilized to compute 
shell displacement and interior pressure responses which are compared' with 
measured values in 1/3-octave frequency bands. Structural damping estimates 
are based on linear viscoelastic theory. Various 1/ 3-octave band averages 
are defined for computing other frequency- dependent parameters for the 
system. Rather good overall agreement between theoretical and experimental 
results for shell response is achieved when the non-ideal characteristics of 
the 1/3-octave filters are accounted for. On the other hand, agreement for 
interior pressure response was somewhat less satisfactory. A detailed 
discussion is given for several possible sources of discrepancy. 
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INTRODUCTION 


Determination of the response of elastic structures to distributed 
random pressures is of fundamental concern in many engineering applica- 
tions. It assumes particular significance in aerospace environments where 
acoustical energy from several sources^ exerts a profound influence on the 
dynamic response of the structure, as well as on its interior components. 
These pressures are usually random in time and are spatially distributed 
in a manner which is dependent on the phase of the launch trajectory. In 
general, the spatial distribution is not random, but is a deterministic 
function which varies with the speed of the vehicle, as well as with its 
position in the trajectory. 

Barnoski, et al. ^ have reported a recent survey of the various methods 
available for the prediction of structural response to random excitation. 
However, most of these methods involve the application of empirical results 
extrapolated from previously available experimental data. The modal, or 
classical method, and the statistical energy method are two- that are more 
generally applicable, and are developed with at least some mathematical 
rigor. Nevertheless, both methods involve numerous simplifying assumptions 
which inevitably limit their use in a given practical case. 

The modal method is generally considered to be useful in low frequency 
regions which include low modal density. Some of its limitations for applica- 
tion to the case of a cylindrical shell excited by random acoustic excitation 
were determined in the initial phase^ of the present program. The results 
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of the remaining work, which are reported herein, represent an investigation 
of the response of the same basic cylindrical configuration over much wider 
frequency ranges which include high modal density and significant acoustic 
radiation. Two different formulations of the statistical energy method are 
employed for prediction and compared with the results of experiments which 
were designed to determine the practical applicability of this method. Thus, 
it is applied to a problem involving an excitation pressure having a space- 
wise deterministic distribution, rather than a diffuse sound field, in order 
to simulate, at least qualitatively, a form of the acoustic fields 
encountered in space vehicle applications. 
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EXPERIMENTAL, ANALYSIS 
Description of Physical System 

Various previous investigations^* have considered the response 
of a cylindrical shell to a spatially diffuse random acoustic excitation. For 
the present case, a non-diffuse excitation is desired. Therefore, the 
physical arrangement depicted in Figure 1 was selected. This is a similar 
system to that which was utilized in our previous effort^, except that more 

elaborate calibration procedures are necessary at higher frequencies, and 

; 

the cylinder was capped in this case. The acoustical speaker was chosen 
to provide a reasonably effective area of excitation, yet small enough to 
minimize the computations required to obtain theoretical numerical 
results. A photograph of part of the actual apparatus is shown in Figure 2. 

The entire apparatus, which was designed to perform several related 
experiments, was composed of six parts: test fixture, excitors, excitation 
sources, transducers, analyzers, and recording devices, all of which are 
shown in the schematic of Figure 3. 

The test fixture was a thin-walled aluminum cylinder whose physical 
properties are given in Table I. This cylinder was bolted (by a welded-on 
ring flange) to a heavy steel plate at its bottom and similarly to a heavy 
3/4-inch aluminum disk plate at its top, which put the cylinder in an e. ficLo ?>■£■&- 
and nearly fixed-end configuration. Wood baffle plates of 1 /4-inch 
thickness separated the interior space from the end plates. An isomode pad 
was used to cushion the steel base plate on crossed steel I-beams, which 
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TABLE I. PROPERTIES OF TEST CYLINDER 

p s = 2.59 X 10" 4 lb-sec 2 /in 4 

h s = 0. 020 in. 

a = 1 2 . 42 in. 

£■ = 30. 0 in. 

f r = 2,640 Hz 

f c = 23,300 Hz 

Material - 6061 -T 6 Aluminum. 



in turn rested on a concrete floor. The function of the isomode pad was to 
eliminate excitation of the cylinder through floor vibrations. A 3/4-inch 
hole was drilled through the center of the aluminum top plate of the test 
fixture to allow access to the interior of the cylinder, so that a study of 
the characteristics of the internal acoustic field was possible. Figure 2 is 
a photograph of the test fixture in the configuration in which the acoustical 
speaker was used as an excitor. 

Two excitors were used. The primary one was ah 8-inch "hi-fi" 
loudspeaker which was mounted in relation to the cylinder as described by 
the coordinate system in Figure 1, where Xq = 15.00 inches. The plane 
defined by the edges of the speaker cone was parallelto and 0. 85 inch from 
the tangent plane to the cylinder at the excitation center (r, 6 , x = a, 0, xq). 
The acoustical speaker was mounted independent of the test fixture, to the 
concrete floor by a steel support as seen in Figure 2. Measurement of the 
properties of the acoustic field generated by the loudspeaker will be 
described under Calibration Procedures. 

The secondary excitor was a small magnet whose pole pieces were 
parallel and close together (approximately 0.25 inch). The coil of this 
magnet could be driven by an AC power source and the pole pieces could be 
placed in close proximity to the cylinder wall (approximately 0. 10 inch). 
The interaction of eddy currents produced locally in the cylinder wall and 
the AC magnetic field effected a remote point excitation of the cylinder at 
twice the frequency at which the magnet was driven. This magnet was 
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used in only a few experiments which were designed to identify acoustic 
modes excited within the interior air cavity. 

The excitation sources were the signal generators used to drive the 
excitors through a 200 -watt Macintosh power amplifier. There were 
essentially four different excitation sources; a sine-wave generator and 
three sources of random noise. The first random source was created by 
taking the output of an Elgenco gaussian random noise generator, equalizing 
by means of LiTV peak-notch filters to compensate for the frequency response 
characteristics of the loudspeaker, and taping the result on an Ampex 
FR-1800L tape recorder for frequencies below 2. 5 kHz. This 60-Hz to 
2. 5 -kHz nominally equalized random noise signal on tape was called our 
wide -band excitation. The second random source was formed by filtering 
this wide -band taped signal through 1/ 3-octave filters, having standard 
center and half-power frequencies as defined in Table II. The taped signal 
was played through these filters one at a time using filters from 100-Hz 
through 2-kHz center frequencies. This was called our 1/3-octave equalized 
signal. The third random source was the Elgenco random noise generator 
filtered directly by the 1/3-octave filters for center frequencies from 100 Hz 
to 5 kHz. This was called our 1/3 -octave non-equalized signal. 

Five transducers were used to measure characteristics of the 
cylinder under excitation. Three of these were Bentley displacement 
detectors located relative to coordinates of Figure 1, while values for these 
probe locations are given in Table III. A fourth transducer was a B&K 
1 /4-inch diameter microphone located at various points inside the cylinder. 
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TABLE XI. 1 / 3 -OCTAVE FILTER CENTER FREQUENCIES 
AND 1/ 2-POWER POINT FREQUENCIES 

Filter No. Center Frequency 1/2 -Power Point Frequencies 


1 

100 

88 

111 

2 

125 

111 

140 

3 

160 

140 

178 

4 

200 

178 

222 

5 

250 

222 

279 

6 

315 

279 

355 

7 

400 

355 

443 

8 

500 

443 

557 

9 

630 

557 

705 

10 

800 

705 

887 

11 

1000 

887 

1112 

12 

1250 

1112 

1405 

13 

1600 

1405 

1778 

14 

2000 

1778 

2220 

15 

2500 

2220 

2790 

16 

3150 

2790 

3530 


17 


4000 


3530 


4440 




TABLE III. 

TRANSDUCER LOCATIONS 



e 

X 

Probe 

1 

81° 

7.55 in, 

Probe 

2 

204° 

12. 25 in 


Probe 3 


318 


24.40 in 
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The fifth transducer was an Endevco accelerometer located on the top plate 
of the test fixture at various points along a radius at 6 =45°. 

Transducer signals were amplified and analyzed by three analog 
methods as shown in Figure 3: (1) The signals were fed directly into a 
Ballantine voltmeter which has an output proportional to the mean square of 
its input, and this output was time-averaged for 5 to 10 seconds as desired. 

(2) The signals were filtered with the 1 /3-octave filters and then the tirhe- 
averaged mean square signal was obtained with the Ballantine meter. 

(3) The signals were fed into a Spectral Dynamics SD-101 tracking filter with 
a 10-Hz bandwidth and then into an SD-109 Power Spectral Density (PSD) 
Analyzer with variable RC. The latter system allowed relatively narrow 
band analysis. The results of all analysis methods w.ere recorded either on 
X-Y recorders or on oscilloscope camera film. 

Calibration Procedures 

In order to define the spatial distribution of the acoustic field, 
calibrations were performed on the speaker prior to its use in the experi- 
ments. The instrumentation setup for the speaker calibration is seen in 
block diagram form (within the dashed-line area) in Figure 3. The B&K 
microphone was mounted in a baffle which simulated a segment of the 
cylindrical shell of the test fixture. The axis of the highly directional 
microphone was coincident with the shell radius, and the entire shell 
segment could be rotated on this radius. Thus, the microphone could 
effectively measure the acoustic field pressure on the cylindrical surface at 
points referenced to corresponding projected points which lay in a plane 



13 

tangent to the segment at its intersection with the speaker cone axis. A 
photograph of the speaker and cylinder segment setup may be seen in 
Reference 3. 

Mapping of the acoustic field spatial distribution was done for the 
center frequencies of the 1/3-octave filters only, and was thus considered 
to be an average over each respective band. The field was found to be 
essentially symmetric with the axis of the speaker cone; thus, only one 
coordinate (R in Figure 1) was necessary to designate a point located in the 
tangent plane, but coincident with a point on the shell at which the sound 
pressure was 'measured with the microphone. 

In order to obtain the field distribution, a cross -spectral density was 
computed between the pressures measured at R = 0 and those for various 
R # 0. These data were found to consist of real (CO) and imaginary (QUAD) 
parts, and were completely deterministic in space as was expected. Data 
were taken relative to CO of CPSD = 1.0 for R = 0. The general empirical 
equation 

Bn 

k^ = exp {-AqR u ) cos (ttR/Pq) (1) 

was found to fit the CO data well for various values of Aq, Bq, and Pq, which 
were dependent on the center frequencies of the 1/3-octave filters. The 
empirical equation 


kj = Dq cos (ttR/Gq) - Bq 


( 2 ) 
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fit the data for the QUAD part where the constants Dq, Eq, and Gq also were 
dependent on the center frequencies of the 1 /3-octave filters. Values of 
these parameters for various center frequencies are given in Appendix A. 
Relative CO and QUAD plots as a function of R are given in Figures 4a and 4b. 

In order to provide a complete absolute calibration of the speaker 
field, the above- relative distribution must be combined with a power spectral 
density measurement at R = 0. For this, the ratio of the time-averaged 
acoustic power in RMS psi squared to the RMS volts squared of signal 
across the speaker terminals was measured for the three different random 
excitation sources at R = 0. For the cases of 1/3 -octave equalized and 
non-equalized excitation sources, the microphone and speaker terminal 
signals were measured directly by the time -averaged mean square Ballantine 
apparatus. For the case of the wide-band excitation source, the microphone 
and speaker terminal signals were analyzed with the 1/3-octave filters and 
the resultant signals measured with the time-averaged mean square 
Ballantine apparatus. A PSD plot of the speaker output in psi^/Hz for the 
constant wide-band input to the speaker terminals (R = 0) is shown in 
Figure 5. This plot was obtained with the use of the Spectral Dynamics 
equipment using a 20-Hz filter, RC = 3 sec. Circled points on the plot at 
the 1 /3-octave filter center frequencies were obtained from the calibration 
constant for the speaker by using the 1/3 -octave equalized excitation source. 
The averaging effect of these wider-band filters compared with the 20-Hz 
filter may be clearly seen. A list of the speaker calibration constants is given 
in Appendix B. 
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Figure 4a. Spatial Distribution Of Acoustic Field ( Real Part ) 
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The Bentley displacement probes were calibrated at locations on the 
cylinder wall of the test fixture with the use of a micrometer head to measure 
static distance to 10"^-inch accuracy. However, this linear range about the 
operating point used during the experiments was more than sufficient to 
insure good dynamic results to 10 - ^-inch accuracy. The axial microphone 
was calibrated with a B&K 124-dB audio pressure standard. The accelerometer 
was calibrated with a Kistler Model 808K/561T quartz vibration calibration 
standard. 

Experiments Performed 

Three major experiments were performed in relation to shell response, 
along with additional supporting experiments. For the appropriate instrumen- 
tation setup and the given switch locations, refer to Figure 3. The major 
experiments involved measurements of shell displacement and interior air 
pressure response for the following conditions: 

(1) Wide-band equalized excitation source--the loudspeaker was used 
as the exciter, and the output of the proximity probes was analyzed 
with the 1 / 3 -octave filters (switch locations A, a, 1-4, I) and the 
Ballantine meter. 

(2) 1/3-octave equalized excitation source — the loudspeaker was 
used as the exciter, and the output of the proximity probes was 
analyzed with the Ballantine meter directly (switch locations 

B, a, 1-4, II). 

(3) l/3-octave non-equalized excitation source--the loudspeaker 
was used as the exciter, and the output of the proximity probes 
was analyzed with the Ballantine meter directly (switch locations 

C, a, 1-4, II). 

Additional supporting experiments were also performed. Top plate 
accelerations were measured at points along a radius for 1/3 -octave 



19 


non- equalized excitation (switch locations C, a, 5, II). A 10 -Hz resolution 
PSD of the outputs of proximity probe No. 2 and the microphone was obtained 
when using random noise through a 10 -Hz tracking filter as an excitation 
source (switch locations D, a, 3-4, III). Damping measurements for the 
tank were obtained by means of sine wave excitation with the speaker and by 
observing the probe outputs on oscilloscope records. Both 1 /2-bandwidth 
and free-decay methods were utilized. Not all components for these experi- 
ments are shown in Figure 3. Finally, some cursory observations of dis- 
placement probe and microphone outputs were made with excitation by means 
of the electromagnetic coil. This method was used to identify low frequency 
air modes as symmetrical. 



THEORETICAL, ANALYSIS 


General Modal Relationships 

Before proceeding to discuss the details of the statistical energy 
method as applied to the vibration of a cylinder in air, first it is necessary 
to recognize the existence of different modal groups over various parts of 
a wide frequency band. Some of the principles set forth by Manning and 
Maidanik 0 will be utilized for this purpose. 

Figure 6 shows a diagram which depicts many of the modes of the 
present cylinder over a wide frequency range. The general relationship 
utilized for calculating these modes is 

v = |pq a^E(n/a) 2 + (m.Q-rr/f ) 2 ] 2 

+ (1 - v|)(m 0 tr/f) 4 /[(n/a) 2 + (m 0 ir/i ) 2 ] 2 } 1/2 (3) 


where 


v = f/f r , Po = (k?/12a2), m = m + 0.2 

& Q 



cf = E/[p s (l - v|)3 


Note, that following Arnold and War burton^, an effective axial wave number 
mg is utilized for the present case of a cylinder with partially fixed ends. 
For convenience, the vertical frequency scale of Figure 6 has been divided 
into the standard 1/3 -octave bands at the 1/ 2-power point frequencies 
previously given in Table II. 
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Figure 6. Modal Diagram For Cylindrical Shell 
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The modes of the cylinder have been separated into three distinct groups 

representing non- radiating (NR), acoustically slow (AS), and acoustically fast (AF) 
modes. Non-radiating and acoustically slow modes are separated by the straight line 


v = (Cg/Cjg )n 


(4) 


while acoustically slow and acoustically fast modes are approximately 
separated by the curve 


n = (c £ /c 0 )v Re |(1 - v 2 ) 1 / 2 - v[l - (v/ v c ) 2 ] 1 / 2 j 1 


(5) 


where 


= fc/f: 


r * 


f c ~ 


D = 


Eh' 


2Tr(D/p s h s ) l/2 5 12(1 - v 2 ) 


It should be recognized that additional, acoustically fast modes occur at 
higher frequencies outside the present range of interest. 

Modal density for the cylinder obviously becomes quite high within 
the frequency range considered in Figure 6. For the sake of information, 
a modal density count for total modal density and density of acoustically fast 
modes are compared with theoretical predictions in Figure 7. Theoretical 
values are based on Eqs. (67) and (68) of Bozich and White-’, which are 
valid for a simply- supported cylinder. There is obvious disagreement for 
the total modal density which may result from the extrapolation of Eq. (3) 
to the case of a cylinder with nearly fixed ends. Therefore, in the present 
work, modal density was based on actual modal count obtained from Figure 6 
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in order to be consistent. Thus, for each 1/3 -octave band, the number of 
each type of modes l2AF» -*-2AS, arL< ^ ^2NR were counted directly from 
Figure 6, along with Eqs. (4) and (5). 

Power Balance Equations 
General Concepts 

The general formulation of the statistical energy analysis 

involves equations for power flow in and out of designated subgroups of a 

2 ft 

system. Barnoski, et al. , and Ungar° have presented recent summaries 
of the fundamentals of the method and are careful to. point out the existence 
of numerous limitations on its use. Further discussion of these limitations 
has been presented by Zeman and Bogdanoff^. For our purpose here, it 
will be convenient to repeat the assumptions given by Reference 2 as 
required in the application of the method: 

(1) Modes of each substructure of interest must be grouped 
into similar sets. 

(2) Coupling between modes in a group is negligible. 

(3) Coupling between groups is conservative. 

(4) Modal damping is light and modal response is mostly 
resonant. 

(5) The power spectrum of force is approximately constant 
over the bandwidth of interest. 

(6) Kinetic energy is evenly divided among modes in a set. 

(7) Kinetic energy in coupling must be small compared to 
modal kinetic energy. 

(8) The coupling factors between modes is constant and not 
strongly frequency dependent near the resonance condition. 
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Use of these assumptions will become apparent in the development 
that follows. 

Detailed Derivation 

One possible form of power balance equations for the vibration 
of a cylinder in air has been presented by Conticelli^ and Bozich and 
White^. However, their equations do not directly contain terms which allow 
for a non-diffuse excitation, or for non-radiating modes. Therefore, we 
will first present a similar set of equations which do account for such 
additions. These equations will be referred to as the separate group theory, 
for reasons which will become obvious. Then, as a result of discrepancies 
which resulted between this theory and measured values in part of the 
frequency range, a second, slightly modified theory is developed. It will 
be referred to as the percentage theory . In particular, the non-interaction 
between AF, AS, and NR structural modes of the shell is considered in an 
alternate manner. Both, however, are slightly different modifications of 
the same basic statistical energy theory. The development will proceed 
from several earlier references on the subject. 

For two normalized linear oscillators^ having instantaneous 
velocities u^ and u^, power balance equations can be written as 


( 6 ) 


Pi( u l) " §12^1 " ®2^ ~ Pl( u l) “ r 12^ u l) “ < u 2> ] ~ ( f l u l) 

P2( u i) “ §21^2 “ 0 l) = ^2 “ r 21 K u |) " ( u l)^ = ( f 2 u 2) 


(7) 
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The first term in Eq. (6) represents the average internal dissipation. 

is a normalized damping coefficient of the first oscillator, and is 
the average energy (per mode) of the first oscillator The 

second term represents the power flow from mode 1 to mode 2; g^ is the 
power flow coefficient; and r-^g is the coupling coefficient between mode 1 
and mode 2. The extreme right-hand side of Eq. (6) represents the work 
done per second, or the input power on mode 1. The terms in Eq, (7) have 
similar meanings. For "gyroscopic" coupling, which holds^ for shell-air 
coupling gj2 = §21* r 12 = r 21> an d t ^ ie dissipation in the system is independent 
of the coupling which is, in this sense, conservative. More general 
coupling equations will be discussed later. 

For our present problem, the above equations are extended to 
groups of modes so that the subscripts 1, 2, and 3 are used for the outside 
air, the cylindrical shell, and the interior air, respectively. However, 
group 2 is further separated into NR, AS, and AF modes. In the separate 
group theory, essentially five groups are then defined, i. e. , 1, 2NR, 2AS, 

2AF, and 3. In this case, we may write Eqs, (6) and (7) accordingly as 
power balance (per unit bandwidth) equations. For the cylinder 
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( E 2AS e 3\ 
^2 E 2AS + UT l23AS n 2AS ^ n 2AS " n 3 j 


E 2AS in 

+ WT| 21AS n 2AS P 


n 2AS 


2 AS 


(9) 


“ r| 2 E 2NR P 2NR 


( 10 ) 


Note that for non-radiating (NR) modes -the acoustic coupling is considered 
negligible. Similarly, for the interior air 


wr l3 E 3 " <or *23AF n 2AF 


/ E 2AF 

V n 2AF 



( K ZAS 

" COT *23AS n 2AS \ n 2AS 


f3 

n 3 




( 11 ) 


A similar equation for coupling to the external air is not required since it is 
assumed to be of infinite extent. These are the governing equations for the 
separate group theory. 

Now, consider a slightly alternate method of arranging the 
governing equations. Consistent with the concepts of the statistical energy 
theory, the modes of each separate medium (1 or 2 or 3) are assumed to 
possess approximately the same energy level. In particular, we shall assume 
that the shell structural modes, AF, AS, or NR are approximately of the 
same energy so the structural coupling between these modes can be neglected. 


We have: 



( 12 ) 



E 2AF _ E 2AS 
n 2AF n 2AS 


E 2NR 

n 2NR 



(13) 


Since the exterior domain of air is assumed to be of nearly infinite extent, 
there is only radiation outward from the shell with negligible reflected 
radiation 


<uf> = 0 (14) 

The power equation of medium 1 thus is decoupled from the shell-interior- 
air power equations. Equations (6) and (7) yield the following average power 
balance equation (per unit band width) for AF, AS, and NR modes, respectively: 


to1 l2 E 2AF + “ r l23AF n 2AF 




+ WT, 21AF n 2AF 


£2 

n 2 


= P' 


IN 

2AF 


(15) 


“’HZ^'ZAS + WT l23AS n 2AS 




+ w1 l21AS n 2AS 


E2 

n 2 



(16) 


WT l2 E ZNR “ P ZNR 


(17) 


Similarly, the average power balance equation for the interior air modes is: 
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■rp / E 2 E 3\ 

“^3 E 3 ‘ a)T l23AF n 2AF ~ —j 


“^AS^AS 


e 2 e 3 


Equations (15) through (18) are the basic governing equations for 
the percentage method. Consider now a slight rearrangement of these 
equations which makes the name of the method more obvious. The total 
energy of the shell is 


E 2 " E 2AF + E 2AS + E 2NR 


The total number of shell resonant modes per unit bandwidth is 


n 2 “ n 2AF + n 2AS + n 2NR 


Only reasonant modes are counted, since except at low frequencies, their 
contribution to energies outweighs other modes provided that they are lightly 
damped. If Eqs. (15) to (17) are summed (with equal weight), one obtains 


/ E 2 E 3 \ / E 2 

^ 2*2 + ^ 23*2 [^'^) + 


- P IN _ pIN , pIN p IN (21) 
7 *2 P 2AF + ^2AS + f 2NR l* 1 1 


where the coupling coefficients are 


n 2AF n 2AS . , 

^2j n , r *2jAF + n , ^jAS ’ j = 1 or 3 
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Equation (22) expresses each part of the radiation as a portion, or percent, 
of the total radiation. This is similar to that presented by Manning and 
Maidanik^ , 

The determination of number of modes I 2 AF’ ^ 2AS » -^2NR 
in each given bandwidth Aw centered at « is described in the previous section. 

For the shell, the respective modal densities are then 

n 2AF = ^Af/^’ etc * (23) 

For the' number of interior air modes, we will employ those for an 
equivalent rectangular room^ (which is only an approximation for the 
present case) 



This expression is valid for a rectangular room at high frequencies, and 
is an approximation based on the more general expression, Eq. (3. 4), on 
p. 86 of Reference 13. Note the difference in notation, and that a differenti- 
ation needs to be performed to get Eq. (24). 

The space -average generalized displacement' for soundwaves in 

* • 11 
air is ■ 



(25) 


' where 4> is the space-averaged velocity potential, while, for small periodic 
disturbances, the space-averaged pressure is: 
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P = ~ Pq -|^ = sin ( <ot )> 4* = $ cos ( ut ) 


( 26 ) 


Thus, 


* 7 / S n3 (w)V 3 

E 3 " Po V 3 <#A> = “ £ ~ 

A 


Po c o 


(27) 


as 


(p^Au 2 S p3 (co)Ao) 


(28) 


Similarly, 


Eg = M { y 2 ) /kco = Mw^Syta) 


(29) 


where y is the space-averaged shell displacement. The space-averaged 
pressure spectral density of the interior, S p3 , and the space -averaged 
shell displacement spectral density, Sy, can be calculated from Eqs. (27) 
and (29) after Eg and E 3 are determined from the power balance equations. 
Input Power 

Previous investigators have usually considered only one dominant 
form of energy loss in impedance relationships which affect input power 
expressions. In the present case, where energy losses vary considerably 
from one part of the frequency range to another, a more accurate approxi- 
mation will be derived. We start from the normalized power Eqs. (9. 14b) 
and (9. 15b) of Lyon and Maidanik^ * : 
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N, 


+ X §mr(®m “ ®r) = Pm^m 


(30) 


N. 


Pr(4 r ) ~ Smr^m " ®r) “ Pr^i 


(31) 


m 


which are the m*'* 1 shell equation and r** 1 interior air equation {effects of 
the exterior air domain will be added later), and 


Pm e m - ( F mS m ) “ ((MM m ) I P<i>m d S * s m) 

A 


(32a) 


where 


F m = IMM m l 4/2 / P<l> m dS 

A 


(32b) 


and (j> m is the m^ 1 shell normal modal function. Since 


<s^> = <s^>; W(X) = £ 


p s h g \l/2 
M 


(q?> = (pVe r /M)(qr Z >; ^(x) = £ q^ r (T) 


m 


M 


f P s h s 

J ^m • M ^n ]y]; ^mn ” e m d mn’ 


m 


A 


= — for simply -supported panel.! 


jNote that this expression will be only approximately correct in the present 
case. 
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/ +A d i = V ^r 6 rk5 e r = ^ if r and k jt 0 
V 


we have 


E s = T— 
5 Ago 


M( " 2> 


(33) 


E a=^P0 v <i 2 > = -t ? P0 Ve r <#> 


(34) 


Now consider a sum of — MX [Eq. (30)] with respect to m. The 

Ago 

first term yields the dissipation, P m E s * The second term is 


M I E Tmr( 0 m “ 0 r) = R rad 


m r 


l < *&> - n s i 

m r 


<q?>' 


N r J 


R rad 

Eg. 

M 

NT 


E 


(35) 


Equation (35) is a generalized definition of rj? as a modification of the 
definition, Eq. (9.19) of Lyon-Maidanik"*- * . The input power in the summed 
equation is then clearly 


N s . . N s 


Pf = I Wm = s £ 


Ago 


m 


(36a) 


m 


where 


F 


m 


1 

Ve m M 



(36b) 


By Fourier transform, similar to p. 32 of Robson^: 
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1 03 

x(t) = — J xfaje^dco 

-co 


( 37 ) 


°o 

x(co) = J* x^Je'^dt 

-oo 


(38) 


Also, the inverse transform of Eq. (8), p. 47, of Robson^, yields 


1 00 

(x(t)y(t + t)) = R xy ( T ) = y J x*ye" lwt dw/2ir 


(39a) 


-oo 


(x(t)y(t)) = R xy (0) =— J x*yd go/2tt (39b) 

-oo 

Consequently, for two narrow bands of width Ago centered at frequencies ± go 
(as complex transforms are used), we have 


(Aco)C xy (oo) = (x(t)y(t)) J x * (go ) y (co ) doo 

-oo 


= (Aoo) 
“ 2 ttT 


.00, Ago 


.go, Aco 


|_S*(Go)y(Go) + x*(-co)y(-w) 


(Ago) 

2ttT 




.ta, Ago 

x*(oo)y(Go) + 


co. Ago 

S(co) y*(co) 


= (Aco)ReS xy (co) 




go, Ago 

x*(co)y(oo) 


(40a) 


Upon division of the last term of Eq. (40a) by Ago, the power spectral density 


results: 
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Sxx(“) = RbS^w) = 


■j, 


(40b) 


while 


_w,Aw 


J_ 2 

2-tt * T 


_w , Aw 


1 - 


Cxy( w ) = E s *(“)y( w )5 = JT -; * 7 pRe[x(w)y*(w)] = C xy (f) (40c) 


2tr 


C xy (f) is the co-spectral density given in Eq. (6.51a), p. 274 of Bendat- 
Piersol^. Similarly, the quadrature spectral density is 


/ i \\ 00 

d *y<“> \ x(t)y K ' zr)) - = 2^ _/ ** y 


dw 


w, Aw 

= - ^ • |jU£*(w)y(w)] (40d) 


The cross-spectral density is, therefore 


■SxyM = \ »(«> y(») = Cxy - i Qxy («e) 

The Fourier transform of Eqs. (9. la) and (9. lb) of Lyon-Maidanik^ 


yields 


N 


r B rm^r 


m 


, T ~ . . > 

r Zjn(w) ^-^(w) 


S ^ — “ ICO s 


m 


(41) 


N" q x? % 

~ _ s B xxn s ui 

ir - E ~pr. , 

m ^r' 05 ' 


= 0 


(42) 
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Z m (w) - + iP m oo, Z r (co) = - (s?* + i(3 r a) 


(43a, b) 


Using Eqs. (36a), (40a), (41), and (42), we have the power input 


P? = ^ 2° < r n y 


m 


, , 7 N s w, Aw 

l (*£An) 

m 


2 it T 


BP p +P gr s 


pIN 


s 1 + a radc 


(44) 


where a rad is a radiation correction factor to be determined later. Employing 
"integral approximation" and residue theorem to evaluate band-average 
values for a lightly damped system analogous to White and Powell^, one has 


-O), Aoj 


pIN,M 


M 


L s 2ir T 


m 


Zm(co) 2 uT 


N P 


m 


F^F 


m-*- m 


(2Aco) 


(45a) 


~ 2 Mm A2ns ^n> s S Po (aj) 


(45b) 


wbiere 


n s = n 2AF’ etc *> ( jn> s = (j^AF’ etc * ’ respectively 

Equation (45b) is the input derived by White - Powell Thus, the present 
problem requires the further evaluation of the radiation correction factor 
a rad> which will be done shortly. 

The joint acceptance in Eq. (45b) can be derived from Eqs. (45a), 


(32b), and (40c) as: 
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< J ^) C =-4- y 


3 a 2 i s 


r /*- pv 

IJ-'J S Po 


x' ) ^ m { x " ) dx 1 dx " 


(46a) 


A’ A 1 


1 Ns r r 

= — E / / EGpip"/Sp o ] Ml’lWs") d x‘ d s" ( 46b ) 

A^I S - m 


A* A* 


p' = p(x',t), p"=p(x”,t), I s = n s Au) 


(47a, b, c) 


The imaginary part of Gpi p iiin Eq. (46b) has no effect due to its anti- 
symmetry with respect to x’ , x". Let 


g p'p 0 =fc a <£ » “ )s p 0 (aj) " jk i ( - 5 ^ s p 0 m 


It can be shown (Appendix C) that 


C p ,pM = ReG plpll = Ek R (x', w)k R (x",w) + k I (x , ,(o)k I (x ,, ,u)] S p (to) (49) 

Consequently, for simply- supported panels and axisymmetric excitation with 
respect to an axis through the center of the speaker and normal to the 
cylindrical shell at (r, 6, x) = (a, 0, xq) 


On). = -At Z <4n + 4n> 


s A^I S -m n 


The first term in the parentheses of Eq. (50) is: 
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r f 0 mx r^mx 

W = J k R (x,o>) cf) mn (x) dx = 4a J J k R (g, 0, u) 

A' 0 0 


cos 


mir| 


j cos (n0) d£ 


muxo 

d0 sin — (51a) 


where 


^mx ®-mx’ ®mx ** s ’- n 


(51b) 


B f 


kp(£,0; co) 2 e -A 0 R cos 


ttR" 


for R < = 8.0 inches 


(51b) 


as given previously by Eq. (1) for the present system. Note that, based on 
experimental data, k R is assumed to be negligible for R > 8. 0 inches where 


R 


= + a^ sin^0 ; £ - 


x - x f 


(51c) 


and Aq, Bq, and Pq are given in Appendix A. Further, 


/ f mx f ^mx 

ki(x, co) 4>m n (x) dx = 4a J J kj(£, 0; co) 

A 0 0 


cos 


( cos ( n 0) d£ d0 sin(^ — (52a) 


where 


kj(£> 0; co) 5 D 0 cos^^Sj - Eq 


- ) 


and Dq, Gq, and E n also are given in Appendix A. 


'0 


(52b) 
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Returning to the correction factor in Eq. (44) with the help of 
Eqs. (41) and (42), we have for N g 0 (all terms are zero when N s = 0). 


^ad 


+ 1 


Re 




. co. > Aco‘ 


I 

m 


„ "r 

Re < y y x 

k. r Z m<“> 


-oo.Aco 


-co, Aco 


Re 


N, 


-co, Aco-i 


Re< 


s ~ ~ 

■j 

Lm 


A mm 


N s 

I i — k#m 

m r 


Z;(co) Z^(co) 




(53) 


Again, from integral approximation, but summing first with respect to r.in 
a procedure entirely analogous to that for Eq. (11-1) of Lyon and Maidanik^, 
the first double sum in Eq. (53) is 


a = Re < y V 
m I A A 


N s N r -B^^s F^Moo' 


• co, Aco 


mr m m 


2" B mr 


R. 


m r Z^(co) Z'^(co) 
n r (co) 


. iyr N r(“) , . r 


r^/ “rn' 
2 


Re 


N 


co, Aco 


z 

m 


s ~ ~ - 

F «t# * 

ih s m 


2Aco 


s m _ ~ lwS m 


it ‘"'rad it co , . , 

2 ’ A coM * 2 ’ A^ (T] 23 + 


(54) 


Considerable effort would be involved in determining a k if k # m, which is 


JN b N r B mr B kr s k F m ( -ico ) 

a k = Re 1 y y - — — 

l k r Z;(co)Z^(co) 


i Aco 


Re 


N s 

I 

m 


.co, Aco 


F” 1 


m s m 


(55) 



40 


N 


s 


a rad = 2 0]*. 
k 


{ 56 ) 


For N s ^ 0, the minimum possible value of a k may be a m while the maximum 
possible value may be N s a m . As an estimate of the correction for later 
discussion, we shall use 


1 + a. 


ad AF £ 1 + 


N s + 1 


^ , . ( l ZAT + 1 

7 a m “ 1 + \ 2 


) 


^23AF + ^ZIAF^ 2 ( 57 > 


Am 


Note that a rac j is negligible for other than acoustic fast modes, 
s 

Coupling Factors 

The radiation coefficient R rac j for a flat panel was derived by 
17 

Maidanik with some simplifying assumptions. However, there is some 

ambiguity and apparent inapplicability. Some revised forms were given 

1 2 

in Crocker and Price which do not seem to be entirely consistent with 
Maidanik's statements,- thus requiring new interpretation. Nevertheless, 
a detailed mathematical rederivation of Eqs. (2-39a, b, c) of Maidanik will not 
be performed here due to limitation of time available and may not necessarily 
be needed. Therefore, we apply the given forrris to the present case directly. 

The coincidence frequency £ c for our experiments is 23, 300 Hz, which 
is much higher than the maximum frequency in our tests. Thus, we shall 
apply equations for this range f < f c , i. e. , Eqs. (2. 39 a, b, c) of Maidanik* 7 . 
This is based on Maidanik's statement above Eq. (2,38) that "We shall 
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therefore confine our attention to frequencies below the coincidence frequency. " 
In Maidanik’s notation list, he defined kp. as panel wave length and fp as 
coincidence frequency of the panel, which will be replaced by f c . 

Equation (2.39a) of Maidanik is then 

■R-rad^g ~ Pg c o {^ a ^^a^-c^ ^ (^r 

for kp > k a , - i, e. , Cp < eg or AS modes (58) 

where 


a = \l f/f c 


(59) 

Gj(f/f c ) = (4/ir4)(l - 2a 2 )/a(l - a 2 ) 1/2 , 

f <-£ r 
~ 2 c 

(60a) 

Gi(f/f c ) = 0, f >~f c 


(60b) 

c ~ J ± k = “ 

P “ k ’ a c 

v K P c o 


(60c) 

G ? (f/f c ) - 1 _ {(1 - a z ) A [(1 + a)/(l - 

2 (2ir) 2 1 

a)] + 2aj - / (1 - a^) 2 ^ 2 

(61) 

X 

K C f 

1 c 


(62a) 

X = — 
f 


(62b) 


The factor 2a 2 in the first term inside the brace of Eq, (58) was first corrected 
by Crocker and Price ° and confirmed by Maidanik through a private 
communication. The range of applicability has been modified and is con- 
sistent with all known usage in the literature, i. e. , for AS modes. 



In our case, the perimeter of an equivalent plate is 


P r = 4ira + 2i (63) 

Equation (2.39c) of Maidanik^ is then 

R rad = A P0 C 0 [1 - £ / £ c l' 1/2=R rad Ai . £or k P < k a. £ -e., 

i. e. , Cp > eg or AF modes (64a) 

Maidanik's branching condition of f > f c is inapplicable since the square root 

I p 

becomes imaginary, while the modification by Crocker and Price is inconsis 
tent with the stated range f < f c . As in our problem £ << f c , Eq, (64a) 
yields 

R rad AF = A Po c O for AF m °des (64b) 

which has been used for panels and is apparently justified. Similarly, 

Eq. (2. 39b) of Maidanik could be for kp = k a which shows a discontinuity in 
radiation coefficient from AF to AS modes. In any event, this equation has 
not been used in practice and will not be discussed further. 

As pointed out by Crocker and Price^ 2 , Eq. (58) applies for AS modes 
which are edge modes, thus proportional to its effective edge. The resistance 
for radiation from the cylindrical shell to the exterior through AS modes is 
then 


_3/4 space 2 1 / Z space 2 I 2ira \ _ 

R 21 AS = R rad AS = 3 R AS =3 [z^TTzl) R rad As 


(65) 



The additional idea of edge modes being inversely proportional to the 

relative value of the space the edge "looks into" was first introduced by 
1 8 

Lyon . It might be justifiable by an integral equation approach, as the 
contribution or residue of a source at the edge is directly proportional to the 
solid angle of the space it looks into. Thus, the local potential would be 
inversely proportional to this angle. This alternative approach might yield 
an improved result or provide a check on Maidanik 1 s * ^ result based on Lyon 
and Maidanik' s ^ approximate formula. Similarly, the resistance for 
radiation from the cylindrical shell into the interior through AS modes is 



d 1/4 space _ 9T3 l/2 space 
R rad A g ” ^ R AS 


.*( 


Zita 


V 2ira + 2 i, 


R 


rad AS 


( 66 ) 


* \ 1 o 

The AF modes, however, as discussed by Crocker and Price 4 are 
surface modes; therefore, they are independent of the total length of the 
edges and the space they look into. The radiation resistance is approximately 
the same for the cylinder and a rectangular plate, that is 


R 21 AF R 23 A F R rad AF 


(67) 


Finally, the coupling factors as indicated by Eqs. (15), (16), and (35) are 


. - R2j s 

^Js Mco 


( 68 ) 


where j - 1,3 and s stands for AF, AS modes, respectively. For non- 
radiation modes, the radiation is negligible. 
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Loss Factors 

Viscoelastic Structural Damping 

The structural damping depends on the physical mechanism of 
the material during loading and unloading^, There are viscoelastic 
damping and hysteresis damping, macroscopically. For an aluminum beam 
with simple-edge conditions, Baker, Woolam, and Young 20 have found that 
the hysteresis damping is negligible compared with viscoelastic damping. 
Their results are in agreement with experiments which account for the 
contribution of air damping. Thus, a similar theory will be applied in the 
present problem. 

Based on viscoelastic theory, the effective modulus of elasticity, 
E v , and Poisson's ratio, v v are 2 ^ 


E v = 


3k u. 


v r“v 


+ T Pv 


v v = 


K ~ j Vv 


2 k + 


3 


(69) 


(70) 


The effective bulk modulus ky. is assumed to be independent of viscoelastic 
effects; thus 


k v = k (71) 

In terms of Lame constants X. and [i, the elastic modulus, Poisson's ratio, 
and bulk modulus are, respectively: 
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E = 


3k^x 


_ X 

v e 2(X + p) 5 


k 



E 

3(1 - 2v) 


(72a, b, c) 


The viscoelastic Lame constants in terms of the Laplace trans- 
form variable s are 


s n b + . 

. . + b „ 


n 

0 

(73) 

‘ 2(s m a m + 

. . . + aQ) 

ii 

i 

co| c\> 

T= 

< 


(74) 


Eor periodic motion proportional to e^°^, 


s = lt») 


(75) 


For standard linear solids, only first derivatives exist; thus 


K'v 


1 + ico (E ’ / E ) 
1 + iuC 1 


(76) 


which becomes p in the limit of zero frequency. From Eqs. (69), (70), 
(71), (72a, b), and (76), one finds 


E v = E 




2v e ) 2(1 + v e ) 


1 + ico(E7E) 
1 + iioC' 


1 


1_ 

+ 3 * 2(1 + v p ) 


}/Wi 

} 


1 + ico(E'/E) 
1 + icoC’ 


(77) 
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f 1 2 

1 + ico(E'/E)l 

.. 1... 1 

/r i 

1,3(1 - 2v e ) 3 

1 + icoC 1 

2(1 + v e )J t 

/ |3(1 - 2v e ) 


1 

+ 3 


1 + ico(E ! /E)" 

1 __.1 

1 + icoC 1 

2(1 + v e )J 


(78) 


In general, the vibration of a viscoelastic shell can be obtained 
from vibration of elastic shell with simple damping term correction. 
However, for nearly incompressible medium v e -> ~ , the limit of Eq. (78) 

Lj 

yields v v = - = v e as expected, and that of Eq. (77) yields 
c* 


E v - 


1 + itofE'/E) 
1 + icoC' 


(79) 


As an approximation, we shall use actual values of v e for v v and E v given 
by Eq. (79). Then, the equation for the vibration of cylindrical shell in terms 
of normal modes W mn becomes 


1 + ia(E'/E) 
1 + icoC 1 


co 


mn^mn + 


9 2 W 


mn 


8t 2 


= 0 


(80) 


9 

Since ico is equivalent to the operator — , this actually means 

9t 


2 w l E ' 2 ^mn 

W -{■ — to — — 

mn mn e mn 


9 2 W 


+ 


mn 


a 3 w 


+ C' 


mn 


= 0 


9t< 


at- 


(81) 


To get an effective damping, let the solution be 


W 


mn 


e (ico - 6)t 


( 82 ) 


Then, the real and imaginary parts of Eq. (81) yield: 
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<*W " a 0 6 “ w 2 + 5 2 + 36C'co 2 - 6 3 = 0 
a.Q<o - 25(0 - C‘co 2 + 36 2 C'co = 0 

where 

E' 2 

a 0 = -E“mn 


(83a) 

(83b) 


(84) 


For a small damping constant 5, the S 2 term in Eqs. (83a, b) probably can 
be neglected, and 


5 = — [ clq - C'(o 2 ] 

Similarly, Eqs. (84) and (85) yield 

co 2 = “mn " S 2 + 2C'co 2 6 = (O 2 ^ - 8 2 + 2C'co 2 6 - (O 2 ^ 

+ 2C'(o 2 8 


(85) 


( 86 ) 


Substituting Eqs. (82) and (86) into (85) and solving for 25, one finds 


P 


mn 


= 25 = 


c' 

) co 2 

• VE 

/ mn 


1 + (C 1 ) 2 co 2 


(87) 


which agrees with Baker-Woolam-Young's results for a beam. Finally, 
the logarithmic decrement is 


^dec 



co — 

mn 


1 + Ccomn 


( 88 ) 



48 


where 


= EL - c\ C = (C') 2 (89) 

E 

Since these constants depend on thickness, known experimental 
are not applicable to our shell of thickness 0, 02 inch. Therefore, 
logarithmic decrements were measured for our shell in air, results of which 
are shown in Figure 8, Matching of the theoretical and experimental 
logarithmic decrements was made at £ = 200 Hz and f = 4000 Hz in order to 
determine Aj and C as 


Ai =5.33X 1 0- 6 /(fff) 

(90a) 

c = 2.28 x nr 7 ) 

(90b) 


The theoretical curve of Figure 8 was then computed by using these values 
in Eq. (88). The above two frequencies were particularly appropriate for 
matching the damping data, as preliminary investigations showed no acoustic 
fast modes in both frequency bands, so that the air damping due to radiation 

r 

is relatively insignificant at these frequencies. Further, the resulting 
theoretical damping allowed good general agreement with experimental 
results for shell response, and plausible explanations for some discrepancies 
to be discussed later. 

Since only resonant modes in the narrow frequency band 
are needed in the prediction, o) mn = co. Thus, the effective damping 
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coefficient for all modes in the frequency band centered at is 
approximately 


Ps - Pz = 2 5 


A i ^ 2 

1 + Cu 2 


1 + (C') 2 to 2 


(91) 


When C ! = 0, Eq. (81) reduces to second order and Eq. (91) yields (E'/E) 

“mn’ w ki c k i s the correct damping coefficient. Comparing Eqs. (7) and (21), 
the loss factor r^ s is then 

*l s = ^2 = (92) 

which is non-dimensional, as it should be. 

Air Damping 

When there is negligible dispersion of sound due to wind and 
turbulence, there are two properties of the medium causing combined attenua- 
tion of a wave propagating in a "free space. " The first is the molecular 
absorption and dispersion in polyatomic gases involving exchange of trans- 
lational and vibrational energy between colliding molecules, and the second 
is dissipation due to viscosity and heat conduction in the medium. The 
attenuation due to the first is and that due to the second is a. 0 , which is 
the sum -of a , due to viscosity, and cu^, due to thermal, conductivity. The 
total theoretical attenuation is 


a A = a M + a C 


(93) 
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Graphs are available to determine aj^ and cxq directly and 
indirectly (see Figures 3d-l and 3d-2 of Reference 22). However, some 
difficulty was encountered, and the resultant damping appeared to be of 
wrong magnitude to check the measured pressure spectral density at 200 Hz. 
Fortunately, ah empirical equation which describes the measured values 
of Knuds en and Harris with good accuracy! for relative humidities above 
30 percent and at temperatures near 20°C is also given, and is valid for 
our experimental conditions.. This equation states 


/ f 0. 085 c „-i 

aA ~\iooo) 20 + 4 > t 


(94) 


where f is the frequency in Hz and 


c^T = <j>20^ ■*" 9.067) AT (non-dimensional) 


(95) 


$20 * s humidity at 20°C and AT the temperature difference from 

20°C. Since Eq. (94) is insensitive to small differences in temperature and 
humidity, we have taken AT = 0 and <j>20 = 0.5; thus, the denominator is 
approximately 20.5. 

-a A x 

The attenuation factor is given as e A ", hence, a one -dimensional 
theory is used to correlate the damping coefficient, i. e. , 


_ e 2 3 P s2 P - 

Cq '3t^ a 9t ° Ox 2 ' 4 9t 2 4 9t 0x2 


( 96 ) 


tThis is true if the unit is 1 / f t instead of dB/ft. 
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The attenuated solution for small damping is 



Therefore 

P a = 2 c 0“A (98) 

Analogous to the relation between Eqs. (91) and (92), the loss 
factor in air is 


r\ 3 = -n a = PJ<* = 2a A c o/“ 


where 


~ { f \ 3/2 0. 085 
aA \1000 Hz/ 20.5 


(99a) 


(99b) 


Based on Figure (3d-3) of Reference 22, the frequency range for the validity 
of Eq. (94) and thus of Eq. (99b), seems to be 10, 000 Hz, beyond which the 
accuracy of approximation would be unknown. In our tests, the frequency of 
interest was below this value, and the use of Eq. (94) or Eq. (99b) seems to be 


well justified. 
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RESULTS AND DISCUSSION 

Comparison of Theory and Experiment 

In effect, results for structural damping have already been presented 
in Figure 8. This, of course, was desirable for describing the method utilized 
to determine material constants. However, at this point, it is appropriate to 
discuss some of its limitations. As can be seen in Figure 8, considerable 
scatter resulted in the experimental data. Such scatter occurs with either 
method used to determine damping. With the free-decay method, scatter 
results from nonuniformity of decay curves which is caused by the beating of 
modes in proximity to each other. On the other hand, with the half-bandwidth 
technique, an erroneous damping usually results from the spatial shifting 
about of modal patterns, as well as coupling between nearby modes. With a 
fixed point for response observation, spatial shifts of modal pattern with fre- 
quency have a strong influence on the damping estimate. A much better 
determination of structural damping in cylinders is badly needed for the pre- 
diction of response to not only random, but to all forms of forced excitation. 

The cylinder and interior air cavity responses, parameters for which 
are given in Figures 9 through 11, are the central r'esults for this entire study. 
Several kinds of data are shown in each figure. Part of the results will be 
discussed in this section, while the rest will be discussed along with additional 
supporting data in the next section. In general, details of theoretical computa- 
tions will first be presented, and then their results will be compared with 
experimental data. 
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Figure 11. Air Cavity Response Ratio 









Figure 9 presents the system energy distribution obtained by means 
of the percentage method only, and is useful for studying some details of the 
statistical energy method. These results are purely theoretical and include 
the assumption of an ideal 1/3-octave filter. On the' other hand, response 
results for cylinder displacement and interior air cavity can more readily 
be compared with measured values. Correlations for these parameters are 
shown in Figures 10 and 11. In these figures, the dashed theoretical curves 
(for both separate group and percentage methods) are based on the assumption 
of an ideal rectangular 1 /3-octave filter, the solid lines represent theoretical 
curves (for percentage method only) which have been corrected for the real 
filter characteristic, while multiple experimental values are given at the 
various frequencies. In Figure 10, measurements are shown for three 
different observation points for broadband equalized excitation through 2500 Hz, 
and 1/3-octave excitation throughout the entire frequency range. In Figure 11, 
the pressure measurements represent maximum and minimum values observed 
along the centerline of the cylinder, as well as the absolute maximum or mini- 
mum value measured anywhere in the tank, for l/3~octave excitation in the 
respective frequency band. 

For convenient reference, the governing equations, on which the 
theoretical calculations are based, are listed in Table IV, along with their 
coefficients and the determining equations for these coefficients. The four 
basic Eqs. (15) through (18) for the percentage method must be solved 
simultaneously to determine the four unknown responses k^AF’ ®2AS’ ^2NR’ 
and in terms of the reference excitation power spectral density Sp^. The 



TABLE IV. EQUATION REFERENCE FOR CALCULATIONS 
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Unknown Variables Determined By Equations 

Percentage Method: (15), (16), (17), (18) 

E 2AF* E 2AS* E 2NR’ E 3 

Separate Group Method: (8), (9), (10), (11) 


Coefficient 


^2 

(92), 

(91), 

(90a), (90b) 

^3 

(99a). 

, (98), 

, (99b) 

r? 23AF 

(68), 

(67), 

(64b) 

r? 21AF 

(68), 

(67), 

(64a) 

r? 23AS 

(68), 

(66), 

(58) 

n 

2 IAS 

00 

(65), 

(58) 

E 2 

(19) 



n 2 

(20) 



n 3 

(24) 



n 2AF 

(23) 



n 2AS 

(23) 



n 2NR 

(23) 




(44), (45b), (50), (51a), (51b), (51c), 
J (52a), (52b), Fig. 6, and 
(57), (68), (67), (64a) 


IN 

P 2AF 

IN 

P 2AS 

P IN 

P 2NR 
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results are then plotted as in Figure 9. Subsequent use of Eqs. (27) and (29) 
then allows a determination of the response ratios as they appear in Figures 10 
and 11. Similar calculations are performed for the separate group method, 
except that Eqs. (8) through (11) are used as the governing equations instead. 

A matrix inversion scheme on a digital computer was used for the 
above solution. Several comments must be made to clarify some of the details 
of this procedure. Calculation of input power coefficients requires the evalua- 
tion of joint acceptance functions by means of Eqs. (51a) and (52a). A numerical 
integration scheme employing a mesh size of 0.4 in. was used for this purpose. 
Further, identification of specific wave numbers, m and n, for each of the 
modes in the respective 1/3-octave bands also had to be included. This 
identification, along with a mode count, was done from the modal diagram in 
Figure 6. All computations then included several parameters which were 
based on some form of 1/ 3-octave average value, which was usually a value 
determined at the band center frequency. These were: 

(1) Spatial distribution of excitation pressure 

(2) Viscoelastic structural damping 

(3) Air damping coefficients 

(4) Radiation coupling factors 

(5) Modal density of interior air cavity. 

It is interesting to look at the results of Figure 9 with one major 
assumption of the percentage method in mind. In particular, the respective 
energies are combined with the modal densities to provide the tabulated 
energy densities per mode given in Table V. It is apparent that the assumptions 
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TABLE V. ENERGY DENSITY DISTRIBUTION 


f, Hz 

E 2AF/ I 2AF 

E 2AS /][ 2AS 

E 2NR /X 2NR 

E 2 /I 2 

Eg / Ng 

125 

-- 

__ 

90. 5 

90. 5 

0 

160 

-- 

-- 

54. 0 

54.0 

0 

200 

-- 

-- 

40.9 

40. 9 

0 

250 


-- 

12. 2 

12. 2 

0 

315 

__ 

-- 

5. 29 X 10" 2 

5. 29 X 10“ 2 

0 

400 

-- 

-- 

3. 95 

3. 95 

0 

500 



0.1192 

0. 1192 

0 

630 


28.3 

0.286 

1. 22 

2. 25 X 10" 2 

800 



0. 138 

0. 138 

0 

1000 

2. 92 

1. 96 

0. 297 

0. 1225 

8. 04 X 10 -2 

1250 

2. 15 

0. 1475 

4. 76 X 10“ 3 

6. 06 X 10" 2 

3. 35 X 10 -2 

1600 

0. 229 

0. 1813 

1.955 X 10" 4 

1. 820 X 10“ 2 

8. 26 X 10-3 

2000 

0. 281 

0.296 

1.812 X 10“ 2 

3. 77 X 10" 2 

1. 635 X 10" 2 

2500 

0. 2156 

0. 1111 

2.46 X 10-3 

6. 03 X 10" 2 

1. 890 X 10" 2 

3150 

-- 

2. 35 X 10" 4 

2. 56 X 10" 5 

1. 023 X 10 -4 

5.41 X 10" 7 

4000 

_ _ 

3. 30 X 10 -3 

3. 28 X 10" 4 

6. 55 X 10" 3 

1. 39 X 10" 6 


5000 
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of Eq. (12) are not very well satisfied by the final results. Nevertheless, 
as will be shown hereafter, the overall comparison of theoretical and experi- 
mental results for the percentage method is still quite good, and is better 
than that for the separate group method. This apparent contradiction simply 
emphasizes the dire need of further study in the application of the statistical 
energy method. 

It must be borne in mind that the theoretical results represent space- 
averaged values. Measured results in Figure 10, which are taken at three 
rather arbitrarily selected points on the cylinder, indicate that the shell 
displacement response becomes reasonably uniform in space only for fre- 
quencies above 300 Hz. This, then, is the most practical frequency region 
where use of the statistical energy method becomes appropriate for the 
cylinder. That is, the actual response values at any point on the tank do not 
deviate appreciably from the average values, so that predicted values are of 
practical use. On the other hand, uniformity of air pressure in th'e-three- 
dimensional interior air space never approaches that experienced by the 
cylinder, and only settles down to about a 6-dB spread above 2000 Hz. It 
appeared in this case that a measurement of maximum and minimum pressures 
was more useful than some average value since such an extreme spread is 
experienced throughout most of the frequency range. Consequently, theoretical 
values which represent a spacewise average should fall within the spread at 
all points. This does not occur near 2500 Hz. Thus, each part of the cylin- 
der and air system independently reaches a state of relatively uniform spatial 
response, which corresponds to a diffuse sound field. This result is probably 
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influenced to some extent by the non- diffuse spatial distribution of the 
excitation. 

The results from the separate group and percentage method (both 
based on an ideal filter) can readily be compared in Figures 10 and 11. For 
the shell response, the two methods give identical results except between 
about 630 and 2500 Hz. This is to be expected since no radiation occurs out- 
side this range. However, within this range, the percentage method can be 
seen to provide a better comparison with experimental results. Similar 
behavior is apparent in the interior pressure response. Except for the single 
point at 630 Hz, the interior air pressure should be negligible below 1000 Hz 
since only non- radiating modes occur below this frequency. However, sig- 
nificant pressures with wide scatter obviously do occur. An attempt at explain- 
ing these observed pressures by means of so-called mass laws proved to be 
futile. Therefore, additional investigations were performed to try to explain 
this discrepancy, as well as the very obvious discrepancies between measured 
values and those predicted by means of theory based on the assumption of an 
ideal, rectangular 1/3-octave filter. The results are given in the next section. 
Supporting Results 

In order to provide a more detailed idea of the frequency variations in 
both shell displacement and interior air pressure, more resolved spectra 
were measured. The results appear in Figures 12 and 13, shown on a relative 
scale. These data are useful to provide a better idea of the respective shell 
and air modal densities. With the apparently wide variations in parts of the 
frequency range, it becomes obvious that measurements with real 1 /3-octave 
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Figure 13. Highly Resolved Air Cavity Response 
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filters that have some finite roll-off can be in error. Further, when the exci- 
tation is filtered, through such a device, low-level noise in an adjacent band can 
often cause a greater response than that in a particular band under concern, 
because of the non-flatness of the acoustical speaker. In order to allow for 
the non-ideal filter characteristic, a filter factor a F and a speaker factor 
a SF now derived. 

The effective input is obtained from a noise generator passing the 
actual filter and speaker. Thus, there is additional power outside the 1/3- 
octave band. It is seen from Eq. (44) that is proportional to the total 
joint acceptance and the integrand of the joint acceptance is implicity related 
to <p(x',t)p(x",t)> which, in turn, is related to p NG from the noise generator. 
We have, in general, analogous to Eq. (40a) 

l r T 

<p(x l ,t)p(x",t)) = — j p(x',t)p(x",t) dt 

0 

oo 

= 2ttT / dt0 

-00 

1 00 

= ■^T J o- F (co,co 0 ) p NG (to,x‘)p NG (6>,x") dw (100) 

-oo 


1 /2 

Therefore, the effective p(co,x) is a^, (o>, ojq)p ng (co, x) . Analogous to Eqs. (44) 

(46a), and (49), we now have 
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Subscript b is for bandwidth number b, and the ^gp 1 are given in Appendix B. 

It is, in our case, sufficient to include the upper and lower neighboring band 
+ and However, the average damping associated with (<) m ^ in bandwidth b 
takes the value at co^. Therefore, it is more accurate to separate the equa- 
tions for +, 0, and - bands and calculate the corresponding mean square shell 
response 

oo 

(y ) - J Sy(“) - Sy(co 0 , co 0 )Aco + ap(co + , w 0 )Sy(cL> + , w + )Aco + 

0 


-H ap(to_, «Q)Sy(a>_, co_)Aco 


(103) 
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which is the Sy(Q) /Sp^(co) using a real filter with factor a F and speaker with 
factor agp. For 1 /3-octave filters utilized in the present study, ap"(co ± , 0) = 

0, 16. Eq. (104) then yields the solid line shown in Figure 10 which repre- 
sents a corrected percentage method, and supports the reliability of measured 
data, as well as the applicability of the theory. A similar correction was 
applied in Figure 11. Additional measurements showed that the significant 
interior pressure response below 1000 Hz was being excited principally by 
motion of the top plate. Acceleration distribution of the top plate at various 
frequencies is shown in Figure 14. The plate apparently is excited by off- 
resonance longitudinal motion of the cylinder in its symmetric modes. 

Further, measurements also indicated that the pressure response was axi- 
symmetric throughout the low range, which tended to confirm its excitation 
by the top plate since axisymmetric cylinder modes occurred in this frequency 
range. Some approximations will now be given to confirm further this 
behavior. 

For a flat top plate moving with a constant or average amplitude W av , 
the solution for its j^ 1 axisymmetric mode is 
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Figure 14a. Top Plate Response Ratio { 200 Hz to 630 Hz ) 








where Xj is the root of J j(X) = 0. Replace Jq^^Wq by for j = 0, 

X. = 0. By averaging over the frequency band, the contribution of each mode 

J 

at the center of the top is 


P 0 co 2 Co(2tt) 
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If we assume W(r, co) is proportional to \ / S-^-, we can then calculate S-yy 

from S-yy as given in Figure 14, and by taking into account the circular domain 

of the top plate, we find, for example, at 200 Hz, with gg = ^3 - 6. 49 X 10" 4 , 

that S w Is = 2, 114 X 10~ 4 . Further, in this range, there is only one 
vv av/ Po 

resonant mode with Xj = 0, the natural frequency of which is less than co'£. 
Therefore, we have 


P3 mx - 


= 4. 64 X 10" 4 = 6. 7 dB Re 10" 4 


(107) 


P0 


This value seems to agree reasonably well with measured data (Figure 11 at 
200 Hz), although it is not exactly the maximum value. 

As frequency increases, more plate modes will appear which could 
excite Xj 0 modes. Strictly speaking, coupling coefficients should be cal- 
culated which become increasingly more complex at higher fr ^encit Due 
to the rapid decrease in plate spectral density, an estimated Sp^ v 
order of 10“® at 4000 Hz, which is much too low. This estimate was 


of the 
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determined by using S-yy( r, «) r for S in Eq. (106a) to calculate S for 

each mode, and then multiplying by an estimated number of fifty axisym- 
metric modes. However, using the actual filter factors, the results at 3150 
and 4000 are in reasonably good agreement with measured data. 

We now consider the possibility of the rectangular room approxima- 
tion as a possible source of error in the predicted values. At center fre- 
quencies of 1000, 1250, 1600, 2000, and 2500 Hz, the value of n 2AS^ n 2AF 
is 2, 1, 2, 3.77, and 5.48, respectively, while no AF mode is present in 
other bands in our present range of interest. Since tjjyr is much greater than 
' ’l AS ( e * £• » ^AS^AF near order of 0. 002 in some cases), the effective 
radiation coefficient when ^AF ^ ® * s a PP rox ^ matel y 


n 


2AF 


n. 


2AF 


,? 23 " n T? 23AF “ n*~~ ^AF 

2 -2 


(108) 


For the above frequencies, Eq. (108) yields 


e 3 = 


£*>0aF E 2 


corj AF E 2 


wJ7 3 + u7? AF ^ Wr? 3 + W77 AF ^ 


(109) 


Since r ? 3 > I 2 are P r °k a kly more "reliable and E 2 yields correct shell 

response density, the doubtful quantity would be N 2 based on the number of 
modes present in the band in an equivalent rectangular room of the same 
volume. If the percentage method or the radiative power flow term is correct, 
the correct value of N 2 should give a good value of E 2 , and the correction 
factor on E^ and thus on Sp^ is: 
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h 

coJ? 3 + wr? AF — 

a N 3 - if < 1I0 > 

“ >r> 3 + wr? AF i 3 

where N3 = and I3 is the accurate number of modes in the band. With 

roots of J 1 (X •) = 0 available in Reference 23 for n — *• m = 0 to 8 and i -» n = 1 
to 5, it is possible to find the exact value of I3 at 1000, 1250, and 1600 Hz. 
For these frequencies. 


I3 = 13, 24, 56, respectively (111a) 

N3 = 17. 9» 36. 5, 75.8, respectively (111b) 

However, with these values used in (110), the correction factor was less than 
1 dB and insignificant. At 2000 and 2500 Hz, the effect is probably also 
small. Thus, the rectangular room estimate for interior aiy modal density 
appears to be a good approximation. 

It next appeared reasonable to question the validity of the radiation 
coefficients as used in the present "percentage" formulation of the statistical 
energy method. It may be that the reverberation effect in the cross- coupling 
term is given in Eq. (36) may be inadequate for this purpose. As a check on 
this possible source of error, we return to Eqs. (30) and (31). Summing 
Eq. (30) with respect to m with weight M/Aco yields 


M ~ n , -2 . M ^ 

Aco £ Pm( s m> Aco ^ 

m m 




- ^ v 8 9 = /F s \ -M_ = 

Aco £ mm \ F m m) Aco P 
m 


(112) 
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The accurate expression of the radiative term is [c.f. , Eqs. (9- 12)> (9. 13), 
(9. 9a), (9.9b) of Eyon-Maidanik^ ^ 
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(113b) 
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It is noted that the radiative coefficient as given by Eq. (9. 31) of 
Eyon- Maidanik * is 

N s (co) N r (co) 
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M 


Z Z ®mr^r’ u m^ 


rad N,,(w) ^ °mr' r' ~m' 

b m r 


M Ns * — N =» 

No 2 ran 
s m 


n r (co) ; n r (co) = n 3 (a>) 


(114a) 


(114b) 


= — 5 — N r (co)&m 

_ TD Ci 


= M ~ B 

2 mn 


n r (co) 


(114c) 


_ p 0 c 0 k a r /■ 1 — — N r (o>) 

= (2 .£ e „| / U + r (i 1 )+ r (5 2 ) d 2; 


s R ' A! 'A 2 


which checks Eq. (11. 1) of the same reference (see Appendix D), if one 
notes that 


N, 


|^ r (5l)^ 2 ) i,r = [^ r (ll) l l J rt2) Nr V^x = i- [sin(k aUl ~^| )/(k a|*l " *z\ )] 
^ -*-1 -*-1 ° 


T — » X 
-2 -*2 


(115) 
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Using this type of approximate processes, it may be possible to reduce 
Eq. (113b) to a more familiar but improved form. This remains to be a 

i i “2i 

future task of research. Lyon-Maidanik has approximated 0 m - Q T by (s y^ ) ; 
however, they have considered only modes of one type. We have changed it 
to - (q r ) in Eq. (3 5)., Both approximations are subject to doubt in a 

general case. 
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RECOMMENDATIONS FOR FURTHER RESEARCH 


It would appear the general results of this study indicate that use of 
the statistical energy method is reasonably adequate for prediction of 
response of coupled cylinder-air systems to excitation by non-reverberant 
acoustic pressure. Prediction of cylinder response is quite acceptable, 
while that for the interior air is somewhat less acceptable. It has not been 
established whether the discrepancies result from the general qopcepts of 
the statistical energy method itself, or from erroneous estimates of some 
coefficient which is used in the method. Thus, considerable additional effort 
is required before a good overall understanding of the capabilities of tfye 
method will be achieved. We will list briefly some of the; steps which may 
be useful to follow in this future effort. 

(1) In order to avoid estimating the errors involved in the use of 
real medium ^-band filters, it would be better to avoid the use of 
typical, commercially -available, part-octave filters. As a 
better approximation, measurements should be made with 
constant bandwidth, narrow band filters, and the results integrated 
electronically over the desired part-octave band. A much better 
approximation of a rectangular filter should result. Unfortunately, 
this was realized too late in the present program. 

(2) Better prediction and measurements of structural damping in 
cylinders is essential. 

(3) The statistical energy method should be applied to studies of 
various kinds of structural elements. It is surmised that some 
of its limitations for use on one kind of element may not be so 
severe on others. 

(4) A further review and more general summary of the statistical 
energy method should be provided. Some recent attempts at this 
have already been set forth. They are particularly desirable 
since much of the earlier literature on the method is usually 
quite terse, and it is fraught with typographical errors and unclear 
mathematical steps. 
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(5) As pointed out at the end of the last section, a more complete 
analysis of the percentage method should be cpnducted. In 
particular, the failure of the detailed results for energy density to 
satisfy the original assumptions, and yet the simultaneous rather 
good agreement of the overall results should be explained. Other 
forms of coupling between various modal groups may be developed. 

(6) The effect of mesh size on input power coefficients should be 
determined. In this summary, all integrations were performed 
with a mesh size of 0.4 in., except for one case. Results for 
Figures 10 and 11 were also computed for a mesh size of 0.2 in. 
at a frequency of 3150 only. It was found that some individual 
joint acceptance integrals changed by as much as 40%, although 
the net effect was negligible on the total for the band. 

(7) In the present experiments, interior air pressures were mea- 
sured at only a few selected points. Measurement at many points 
would allow the calculation of a space average which could more 
directly be compared with predicted results. Howeve.r, in the 
presence of wide spatial variations, the practical use of the 
results is still in question. It would appear that a prediction of 
maximum response would be of more usethan that for the average 
response. 
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APPENDIX A 

SPATIAL DISTRIBUTION PARAMETERS FOR 
ACOUSTICAL FIELD 



SPK. SPATIAL DISTRIBUTION EMPIRICAL EQUATION CONSTANTS 


1/3 -Octave 
Center Frequency 

A 0 

B 0 

P 0 

D 0 X 10 2 

E 0 X 10 2 

Go 

100 

3. 47 X 10 “ 2 

2 . 682 

150. 0 

0 

0 

0 

125 

3. 47 X 10 ~ 2 

2 . 682 

150. 0 

0 

0 

0 

160 

3. 47 X 10 " 2 

2 . 682 

150. 0 

0 

0 

0 

200 

3.47 X 10 ~ 2 

2 . 682 

150. 0 

0 

0 

0 

250 

3.47 X 10 ~ 2 

2 . 682 

150. 0 

0 

0 

0 

315 

3.47 X 10 ' 2 

2.682 

150. 0 

0 

0 

0 

400 

3.47 X 1 0 " 2 

2.682 

150. 0 

0 

0 

0 

500 

3.47 X 10 ~ 2 

2 . 682 

150. 0 

0 

0 

0 

630 

3.47 X 10 " 2 

2 . 682 

150. 0 

0 

0 

0 

800 

3.47 X 10 " 2 

2 . 682 

150. 0 

8 . 0 

8 . 0 

6 ; 0 

1000 

0. 00856 

2.430 

150. 0 

6 . 0 

6 . 0 

4. 0 

1250 

0. 00172 

1.906 

11.0 

15. 0 

15. 0 

4.0 

1600 

0. 00520 

2 . 130 

9.0 

- 10 . 0 

-10. 0 

4. 0 

2000 

0. 3752 

0. 813 

3.0 

- 30. 0 

- 20 . 0 

3.3 

2500 

0. 1685 

0 . 799 

4.05 

60. 0 

40. 0 

2.5 

3150 

0. 0836 

2.330 

5.25 

12 . 0 

4. 0 

2. 5 

4000 

0 . 3092 

0.724 

2.00 

30.0 

- 8 . 0 

1.5 


CO 

/itR\ 

kj^ = exp {-AqR u ) cos j f° r R <8. Ot 

Quad 

kj = Dq cos - Eq for R. <. 8. 0$ 


f For f = 2500 k x = 1 . 00 for R < 1 . 5 
JF or f = 4000 ky = 0. 00 for R < 1 . 0 
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APPENDIX B 

SPEAKER CALIBRATION FACTORS 
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SPEAKER CALIBRATION CONSTANTS 


1 / 3-0 Center 
Freq. 

1/3-0 Band 
into Spk. No. Eq. 
(psi/V cms )2 X 106 

1/3-0 Center 
Freq. 

1/3-0 Band 
into Spk. No. Eq. 
(psi/V cms)^ X 10^ 

100 

0.989 

630 

3.04 

125 

1. 113 

800 

3.98 

160 

1.043 

1000 

5. 66 

200 

0.937 

1250 

1.855 

250 

0.967 

1600 

1.633 

315 

1.026 

2000 

1.056 

400 

1.336 

2500 

0.432 

500 

2.02 

3150. 

0.244 



4000 

0.311 


5000 


0.476 
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APPENDIX C 
DERIVATION OF EQ. (49) 
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DERIVATION OFEQ. (49) 
For narrow-band sampling, consider 


R I 

Pl = Pi COS (cot) + Pj sin (cot) 

(C-l) 

P 2 = P^- cos («t) + p* sin (cot) 

(C-2) 

Define p° as p with 90° shift* 5 . yf e ^ ave 


Pj ~ ” Pj 1 sin (“*) + Pj cos (cut), j = I, 2 

(C-3) 

The time averages yield 


? Pl P 2 = = I <pM ■ p l p 2> = - <P 2 Pi> = Q p 2Pl 

(C-4) 

^ P 1 P 2 = = 2 (p ? p 2" + p l p 9 = ^ P 2 P 1^ = C p 2 Pl 

(C-5) 

Let 


p = p(x,t), p' = p(x',t), p"= p(x",t), Pq = p(0,t) 

(C-6a, b, c) 

By approximating the measured data, we have 


A 

C P Q P = k R(S>“) S p 0 (“)» Qp 0 p •= kx(x,oj)Sp o (td) 

(C-7a,b) 

where 


Sp 0 {«) = Cp 0 p Q = Gp oP() = ^ [(pB) 2 + (p*) 2 ] 

(C-8) 
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For convenience, define p^- = p R ; then, we have 

^P'P 0 = J ^ P R P '0 _ p l p 0 l) = " k I^S >CJ ) S p 0 ( w ) 

^p"p Q ~ 2 ^ P R P 0 ” P I P ?> = “ k j(2i"»“) S p^(“) 

Multiplication of Eqs. (C-8) and (C-9) gives 

k I (x , , Cl ))k I (x ,, co)s|^(«) = \ ^Pr(Po ) 2 + PiPilPo - ) 2 

- <P£p£ + P^P^P^P^ 

Similarly, one obtains 

k R (x , ,co)k R (x",o))S^ (oj) = | 1 p r Pr(Pq") 2 + PiPi(Pfj) 2 

+ + PkPl’PoPo 1 

Adding of Eqs. (C-ll) and (C-12) yields 

k *■ P R P R + P i P I ] ««> + ( pI o ) 2 ^ = Ife I (s' > u)k I {x M ,u) 

+ k R^ , ' w ^ k R^ M »“^ s p Q («) 

Dividing Eq. (C - 1 3 ) by S p ^ with Eqs. (C-8) and (C-5), one finds 

C p i p n = Ik I (x,co)k I (x',co) + k R (x,co)k R (x",<o)] Sp^co), (Q.E.D) 


(C-9) 

(C-10) 

(C-ll) 
(C-12) ■ 
(C-13) 


(C-14) 
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APPENDIX D 

DERIVATION OF EQS. (114b) AND (105) 



DERIVATION OFEQS. {114b) AND (105) 


"We have Eq. {9. 11) of Lyon-Maidanik'*' * that 
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B Lr«W>l + f>r“m> 


mr 


- “r ) 2 + (Pm + P r >(Pxn«? + P r 4 i) 


(D-l) 


By changing summation to integration by assuming there is a large number 
of air modes as done by Lyon-Maidanik and by residue theorem after taking 
advantage of the 8 -function type integrand, we found that 
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(D-2) 


which justifies Eq, (114b) from Eq. (114a). Next, we shall derive Eq. (105). 
From Eq. (9. 6) of Lyon-Maidanik 
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5Z-N r («) 

"“mn 


PQ c O 
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MVe r c m > 


;/ 

Al An 


N r { w ) 


<W£1> W52> dx l d 32 
(D-4) 


Since 1® a two-dimensional vector function on the shell surface, we only- 

have two components of the surface vector; therefore, in place of Eq. (7. 13) 
of Lyon-Maidanik, we should use 


h<5i 


1 


>*r<52>=T6 2 

h= 1 


2 

I 

j = l 


ik.jx^ 


J + e 


- ik j x hj 


(D-5) 


Now, instead of 8 positions in space, we only have 4 positions on the surface 
that the exponent has minimum value proportional to the difference of 

and x 2 such as exp [k x (x^j - x 2 ^) + ky(xj 2 - x 22 )1. Now, we should average 

* ? 7 1 / Z 

with respect to k x and ky over the circle of radius k a = (k x + ky) 

2 2l/2 

= (k a - k^) 1 and average over k z from zero to k a . Therefore, with 6'as 
the angle between the vector k a and xj - x 2 , we have 


+ r ( x l)+r(5 2 ) 
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d0 dkz 

P-6) 


= 1 h-f a J 0 < 1—1 - Szl - k*) dks (D-7) 
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1 

sin (|Xj -x 2 | 

kz) 

k 
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.1 

sin (k a |x 2 - x 2 | ) 
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~ 4 
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Equations (D-7) and (D-8) were derived with well-known formulas given 
by McLachlan (2nd edition) 24 .. If we follow through the factors in Lyon- 
Maidanik’s paper, we can show that 






sin (k| x 


¥l - I ) 


8 

To -X, 


k 3V - 32 


(jl* -i?2 ^Lyon-Maidanik 
(D-9) 


as there were 8 positions in space divided by 64. The factor — should be 

8 

present in Eq. (7. 16) of Lyon-Maidanik and Eq. (2. 18) of Maidanik. 

Equations (D-8) and (D-9) verify Eq. (115). The final result of Lyon- 
Maidanik was correct as they had a factor 2 too large in using their 
Eq. (9. 17) instead of Eq. (D-2) and a factor 2 too small by using Eq. (D-9) 
instead of the correct average Eq. (D-8). It may be stated that Eq. (114c) 
yields Eqs, (2. 15) and (2.21) of Maidanik, which rightfully yields the 
approximate radiation formulas, Eqs. (2.24), (2. 25) in his paper with his 
approximation Eq. (A- 1 ); namely, 


R 


rad 


= A Po c o [1 - 


<V fe a> 21 


• 1/2 


= A 


PQ c 0 


[1 - 


7-1/2 

< V k a } 1 


= A p 0 Cq [ 1 - (kp/k a ) 2 ] “I/ 2 for kp < k a , AF modes 

R rad= 2<Ap 0 Co/3ir i/2 )(ik2/k x ) l/2 

- 2(A Pq c o / 3irl / ^)(hkp/ky)l/^ 
- ( A Pq c o / 3 Trl ' / 2 )t(- ek p / k x ) 1 ^ 2 + (hkp/ky) 1 ^ 2 ] k p = k a 


(D- 10) 


P-11) 



Due to time limitation, however, Eqs. (2. 39a, b, c) of Maidanik have not 
been rederived. However, at least Eq. (2. 39a) was checked and corrected 
by Crocker-Price and is used for AS modes (kp > k a ) in this paper. 



